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We report the superconductivity in Ca-intercalated bilayer graphene CeCaCe, the thinnest limit 
of Ca graphite intercalation compound. We performed in situ electrical transport measurements 
on pristine bilayer graphene, CeLiCe and CeCaCe fabricated on SiC substrate under zero and non¬ 
zero magnetic field. While both bilayer graphene and CeLiCe show non-superconducting behavior, 

CeCaCe exhibits the superconductivity with transition temperature (Tc) of 4.0 K. The observed 
Tc in CeCaCe and the absence of superconductivity in CeLiCe show a good agreement with the 
theoretical prediction, suggesting the importance of a free-electron-like metallic band at the Fermi 
level to drive the superconductivity. 

PACS numbers: 68.65.Pq, 72.80.Vp, 73.22.Pr, 74.70.Wz 


Atomic-layer superconductors (ALSCs), where only 
one or a few atomic layers at the surfaces or interfaces be¬ 
come superconducting [1-8], have attracted considerable 
attentions owing to their two-dimensionallity[5-7] and 
possible existence of symmetry mixed Cooper pairs [9]. 
It is known that in ALSCs fabricated on semiconductor 
substrates, the super conduct ing-transit ion temperature 
(Tc) is in general suppressed in comparison with that of 
bulk [1-5] due to the interference from the substrate ex¬ 
cept for limited cases[7, 8]. In this sense, self-standing 
ALSCs free from the substrate effect have been desired 
to pursue higher Tc. 

Graphene is a single atomic sheet of graphite and now 
a target of intensive theoretical and experimental studies 
because of its various peculiar but attractive properties 
such as the massless nature of carriers [10], the anoma¬ 
lous quantum Hall effect [10], the high carrier mobility to 
drive the ballistic transport [11], etc. Besides these su¬ 
perior properties, graphene has a very important inher¬ 
ent property, namely, self-standing crystal structure ba¬ 
sically isolated from substrate. In this context, graphene 
is regarded as one of ideal materials to realize ALSCs 
free from the substrate effect. Just after the discovery 
of graphene[12], many intensive efforts have been made 
to fabricate superconducting graphene by intercalating 
metals like in graphite intercalation compounds (GICs). 
Ca-intercalated bilayer graphene (CeCaCe) is regarded 
as the most promising candidate because corresponding 
Ca-intercalated GIC (CeCa) has the highest Tc of 11.5 K 
among all GICs[13-16]. While some theoretical[17, 18] 
and spectroscopic [19] studies have proposed a similar 
electronic structure between CeCa[20-22] and CeCaCe, 
suggestive of superconductivity in CeCaCe, no direct ev¬ 
idences for superconductivity, such as superconducting 
gap and/or zero-resistance, have been reported. This 
may be mainly due to difficulty in handling reactive sam¬ 
ples containing Ca, requiring an in situ measurement un¬ 


der ultrahigh vacuum at ultralow temperature [2-5]. 

In this Letter, we report on the observation of the 
zero-resistance state in CeCaCe to directry prove the oc¬ 
currence of supercondutivity by in situ electrical trans¬ 
port measurements. We compare the experimental re¬ 
sults with those for pristine bilayer graphene and CeLiCe 
as well as the theoretical predictions to obtain an insight 
into the superconducting property and mechanism in in¬ 
tercalated bilayer graphene. 

To fabricate intercalated bilayer graphene [Fig. 1(a)], 
we at first prepared a high-quality pristine bilayer 
graphene sheet on a n-type Si rich 6i7-SiC(0001) sin¬ 
gle crystal by heating the crystal up to 1550°C in an 
argon atmosphere[23]. After a short exposure to air, 
the grown bilayer graphene/SiC was transferred into a 
RHEED (reflection-high-energy electron diffraction) vac¬ 
uum chamber where the transport measurements were 
also done. After heating the sample at 400°C for several 
hours in the chamber, we observed a RHEED pattern 
typical of clean bilayer graphene on SiC [Eig. 1 (b)]. Then, 
we deposited Li atoms on the bilayer graphene sheet us¬ 
ing a SAES Getter dispenser. The Li deposition pro¬ 
duced several sharp x ^/3 R30° spots in the RHEED 
pattern [pointed by yellow arrows in Eig.1(c)], indicating 
that Li atoms are intercalated in a regular manner be¬ 
tween two adjacent graphene layers[24] as in bulk CeLi. 
After confirming the growth of high-quality CeLiCe on 
SiC, we then deposited Ca atoms on this CeLiCe sheet 
to replace Li with Ca. During the Ca deposition, we kept 
the substrate at 150°C, slightly above the Li desorption 
temperature of 145°C. We observed that the Ca depo¬ 
sition transformed the \/3 x R30° spots into streaks 
[Eig. 1(d)], suggesting that intercalated Li atoms are re¬ 
placed by Ca atoms [19]. Repeated cycles of Li and Ca 
deposition together with annealing at appropriate tem¬ 
peratures made the ^/3 x \/3 R30° streaks brighter and 
brighter. Einally, we obtained the RHEED pattern in 
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FIG. 1. (Color Online) (a) Crystal structure of CeCaCe 
on SiC substrate, (b) RHEED pattern of pristine bilayer 
graphene. Si(lx 1) spots from the substrate and C(1 x 1) spots 
from graphene are shown. Spots from the buffer layer with 
6a/3 X 6a/3 R30° periodicity are also seen, (c) RHEED pat¬ 
tern of CeLiCe. (d), (e) Same for CeCaCe (d) after the first 
Li-Ca replacing treatment and (e) after several cycles of the 
replacing treatments. Yellow allows indicate a/3 x \/3 R30° 
spots and streaks. 


Fig. 1(e), where we observe sharp a/3 x \/3 R30° streaks, 
indicative of well-ordered CeCaCe- 

After fabricating the sample as described above, we 
then transferred it on a stage for in situ electrical trans¬ 
port measurements with the 4-point-probe (4PP) method 
in the same vacuum chamber[25]. As shown in Fig. 
2(a), the 4PP chip composed of four copper wires in 
100 jam (j) was contacted on the sample, and then cooled 
down to 0.8 K together with the sample. The results of 
temperature-dependent transport measurements on bi¬ 
layer graphene, CeLiCe and CeCaCe are compared in 
Fig. 2(a). The sheet resistance (i^sheet) of pristine bi¬ 
layer graphene shows a metallic character above 20 K, 
but turns into an insulating one at lower temperatures. 
This is usually observed in the transport characteristics 
of graphene and interpreted in terms of the localization 
effect due to the quantum interference [26-28]. The resid¬ 
ual resistance at 0 K is estimated to be 565 ft [Fig. 2(b)]. 
In CeLiCe and CeCaCe, on the other hand, the i^sheet be¬ 
haves as metallic, and their resistance is as low as ca. 10% 
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FIG. 2. (Color Online) (a) Sheet resistance Rsheet of pristine 
bilayer graphene, CeLiCe and CeCaCe as a function of tem¬ 
perature. Inset shows a schematic picture of 4-point-probe 
measurement set-up. (b) Rsheet of pristine bilayer graphene 
in an expanded scale from 0.8 to 60 K. Dashed line shows 
extrapolation toward 0 K. (c) Rsheet of CeCaCe from 0.8 K 
to 20 K, showing = 4.0 K and = 2.0 K. 


of that of pristine bilayer graphene. This is because of the 
increase of the Fermi-surface volume by the career doping 
from Li or Ca atoms and the resultant folding of Brillouin 
zone due to the induced \/3x \/3 R30° superstructure [19]. 
Importantly, the transition of i^sheet to zero-resistance is 
clearly seen at around 2 K in CeCaCe, which is the di¬ 
rect evidence of occurrence of superconductivity. It is 
also noteworthy that CeLiCe does not show supercon¬ 
ductivity down to 0.8 K, and instead exhibits a weak 
localization behavior as evident from a slight upturn in 
resistance at low temperature. Figure 2(c) shows the re¬ 
sistance of CeCaCe at temperatures from 0.8 K to 20 K, 
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FIG. 3. (Color Online) (a) i^sheet of CeCaCe as a function 
of temperature under different magnetic fields, (b) Same as 
(a) as a function of magnetic field for different temperatures. 
Magnetic field was applied perpendicular to the sample sur¬ 
face. (c) (Picture) RHEED pattern of this sample, showing 
weaker \/3 x \/3 R30° streaks than those in Eig. 1(e). Note 
that the sample is different from that in Eig. 1(e). (Plot) The 
intensity of RHEED pattern as a function of horizontal posi¬ 
tion for Eig. 3(c) (upper black dot) and Eig. 1(e) (lower red 
circle). The intensity is normalized by the peak height of cen¬ 
tral (1x1) streak, (d) Temperature dependence of the upper 
critical field /io^c 2 obtained from (a) (squares) and (b) (tri¬ 
angles). 1^0 Hc 2 is defined as the magnetic field where Rsheet 
drops to a half of normal-state resistance. Solid line shows 
the fitting based with the Ginzburg-Landau theory[40]. 


highlighting that i^sheet suddenly drops at around 4 K 
(T°nset _ 4 Q notice in Fig. 2(c) that 

i^sheet starts to dccrcasc even above 4 K. This is prob¬ 
ably due to the superconducting fluctuation inherent to 
low-dimensional superconductors [29]. 

To obtain further evidence for superconductivity, we 
conducted magnetoresistance measurements on another 
CeCaCe sample prepared with a same method as de¬ 
scribed above. In the experiments, the magnetic field 
was applied perpendicular to the surface. Figures 3(a) 
and (b) show i^sheet data as a function of temperatures 
and magnetic fields. As seen in Fig. 3(a), the grad¬ 

ually shifts to lower temperature as the magnetic field is 
increased. A similar behavior of is also seen in Fig. 

3(b). It is noted that the i^sheet curve did not fully drop 
to zero-resistance even at 0.8 K under zero magnetic field, 
suggesting that only a limited part of sample became su¬ 
perconducting and the superconducting paths were not 
fully connected between the current probes in this sam¬ 
ple. In fact, as shown in Fig. 3(c), the x R30° 


streak originating in CeCaCe is weaker in the present 
sample than in that of Fig. 1(e). This in return supports 
that the superconductivity emerges in the well-ordered 
\/3 X \/3 R30° superstructure. To see the temperature- 
dependence of the upper critical field (/io^c 2 ), we plot 
the magnitude of magnetic field at which the i^sheet is a 
half of the normal-state-resistance as a function of tem¬ 
perature in Fig. 3(d). We find that extracted points 
are well aligned linearly as seen in Fig. 3(d), indicat¬ 
ing that these experimental results can be analyzed in 
the framework of the Ginzburg-Landau (GL) theory [30]. 
Numerical fittings with the GL theory show that the in¬ 
plane GL coherence length at zero Kelvin ^(0) is 49 ± 
1 nm. This is comparable to or slightly larger than the 
in-plane ^(0) reported for bulk GeGa (29 - 36 nm)[13- 
16]. This shows a striking contrast with the case of usual 
ALSGs such as In on Si(lll), where the ^(0) ^ 25 nm[4] 
is much shorter than that of bulk (250 - 440 nm)[31]. 
This suggests the strong two-dimensional nature of su¬ 
perconductivity in both GeGaGe and GeGa. 

Now we discuss the origin of superconductivity in 
GeGaGe. While the Tc is lower than that of bulk GeGa 
(11.5 K), the present observation of superconductivity 
in GeGaGe is striking. Because, no superconductivity 
has been observed so far in second-stage GIGs (G 12 A 
where A is an alkali metal) which are a bulk analog of 
intercalated bilayer graphene (GeAGe). A simple esti¬ 
mation indicates that the carrier density in a single car¬ 
bon layer in GeGaGe should be reduced to a half of that 
in bulk GeGa. However, Mazin and Balatsky[17] have 
theoretically proposed the possibility of superconductiv¬ 
ity in GeGaGe despite the reduced carrier density, based 
on the similarity in the topology of Fermi surfaces con¬ 
sisting of a nearly-free electron band at the center of 
Brillouin zone and the tt* band at the zone boundary. 
They showed by calculations that the nearly-free electron 
band is partially occupied in GeGaGe while it is totally 
unoccupied in GeLiGe, indicating the importance of the 
metallic nearly-free electron band for the superconductiv¬ 
ity. This theoretical prediction shows a good agreement 
with the present observation that the superconductivity 
emerges in GeGaGe while not in GeLiGe- Jishi et a/. [18] 
have confirmed the band structure of intercalated bilayer 
graphene proposed by Mazin and Balatsky and further 
predicted the Tc of 4.7 K for GeGaGe- This shows a 
good quantitative agreement with the present observa¬ 
tion of Tc = 4.0 K. This strongly supports the theoreti¬ 
cal proposal[17, 18] that the metallic nearly-free electron 
band at the center of Brillouin zone plays an essential 
role for the superconductivity. One may notice here that 
the calculations [17, 18] were done for ideal free-standing 
intercalated bilayer graphene while in reality the inter¬ 
calated bilayer film is fabricated on SiG substrate, which 
may affect the electronic structure and consequently the 
superconductivity. In fact, a recent scanning tunnel mi¬ 
croscope study reported that a charge density wave is 
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created at low temperatures in CeCaCe on SiC due to the 
commensurate lattice matching [32]. It should be noted 
that despite these interferences from the substrate the 
Tc observed in CeCaCe fabricated on SiC shows almost 
the same value as predicted for a free-standing CeCaCe, 
suggesting that the superconductivity in CeCaCe is very 
robust. 

In conclusion, we performed in situ electrical trans¬ 
port measurements on Ca-intercalated bilayer graphene 
CeCaCe under zero and non-zero magnetic field. We have 
clearly observed that the resistance steeply drops at 4 K 
and reaches zero at 2 K under zero magnetic field, show¬ 
ing the occurrence of superconductivity with T^^set _ 
4.0 K and = 2.0 K. The measurement under mag¬ 
netic field has confirmed the superconductivity origin of 
the observed zero resistance. The observed is in 

good agreement with the theoretical prediction. The ab¬ 
sence of superconductivity in pristine bilayer graphene 
and Li-intercalated bilayer graphene CeLiCe in contrast 
to CeCaCe suggests the importance of free-electron-like 
metallic band at the Fermi level to ignite the supercon¬ 
ductivity. 

We thank K. Suzuki and T. Sato at Tohoku University 
for their useful discussion. This work was supported by 
the JSPS (KAKENHI 15H02105, 25246025, 22246006) 
and the MEXT (Grant-in-Aid for Scientific Research on 
Innovative Areas “Science of Atomic Layers” 25107003 
and “Molecular Architectonics” 25110010). 

* Corresponding author: 

ichinokura@ surface. phy s. s. u-1 okyo .ac.jp 


[1] S. Qin, J. Kim, Q. Niu, and C.-K. Shih, Science 324, 
1314 (2009). 

[2] T. Zhang et aL, Nat. Phys. 6, 104 (2010). 

[3] T. Uchihashi, P. Mishra, M. Aono, and T. Nakayama, 
Phys. Rev Lett. 107 , 207001 (2011). 

[4] M. Yamada, T. Hirahara, and S. Hasegawa, Phys. Rev. 
Lett. 110 , 237001 (2013). 

[5] C. Brun et aL, Nat. Phys. 10, 444 (2014). 

[6] N. Reyren et aL, Science 317 , 1196 (2007). 

[7] H.-M. Zhang et aL, Phys. Rev. Lett. 114, 107003 (2015). 

[8] J.-F. Ge, Z.-L. Liu, C. Liu, C.-L. Gao, D. Qian, Q.-K. 
Xue, Y. Liu, and J.-F. Jia, Nat. Mater. 14, 285 (2015). 

[9] L. P. Gor’kov and E. 1. Rashba, Phys. Rev. Lett. 87, 
037004 (2001). 

[10] K. Novoselov, A. K. Geim, S. Morozov, D. Jiang, M. Kat- 
snelson, 1. Grigorieva, S. Dubonos, and A. Firsov, Nature 
438, 197 (2005). 

[11] K. S. Novoselov, Z. Jiang, Y. Zhang, S. Morozov, 
H. Stormer, U. Zeitler, J. Maan, G. Boebinger, P. Kim, 
and A. Geim, Science 315 , 1379 (2007). 

[12] K. S. Novoselov, A. K. Geim, S. Morozov, D. Jiang, 
Y. Zhang, S. Dubonos, 1. Grigorieva, and A. Firsov, 
Science 306 , 666 (2004). 

[13] T. E. Weller, M. Ellerby, S. S. Saxena, R. P. Smith, and 
N. T. Skipper, Nat. Phys. 1 , 39 (2005). 


[14] N. Emery, C. Herold, M. d’Astuto, V. Garcia, C. Beilin, 
J. F. Mareche, P. Lagrange, and G. Loupias, Phys. Rev. 
Lett. 95, 087003 (2005). 

[15] E. Jobiliong, H. D. Zhou, J. A. Janik, Y.-J. Jo, L. Balicas, 

J. S. Brooks, and C. R. Wiebe, Phys. Rev. B 76, 052511 
(2007). 

[16] R. Xie, D. Rosenmann, A. Rydh, H. Claus, G. Kara- 
petrov, W. Kwok, and U. Welp, Physica C: Supercon¬ 
ductivity 439, 43 (2006). 

[17] 1. Mazin and A. Balatsky, Phil. Mag. Lett. 90, 731 (2010). 

[18] R. Jishi, D. Guzman, and H. Alyahyaei, Adv. Studies 
Theor. Phys. 5, 703 (2011). 

[19] K. Kanetani, K. Sugawara, T. Sato, R. Shimizu, 

K. Iwaya, T. Hitosugi, and T. Takahashi, Proc. Natl. 
Acad. Sci. USA 109, 19610 (2012). 

[20] G. Csanyi, P. Littlewood, A. H. Nevidomskyy, C. J. 
Pickard, and B. Simons, Nat. Phys. 1, 42 (2005). 

[21] K. Sugawara, T. Sato, and T. Takahashi, Nat. Phys. 5, 
40 (2009). 

[22] S.-L. Yang, J. Sobota, C. Howard, C. Pickard, 
M. Hashimoto, D. Lu, S.-K. Mo, P. Kirchmann, and 
Z.-X. Shen, Nat. Commun. 5 (2014). 

[23] K. Sugawara, T. Sato, K. Kanetani, and T. Takahashi, 
J. Phys. Soc. Jpn. 80, 024705 (2011). 

[24] K. Sugawara, K. Kanetani, T. Sato, and T. Takahashi, 
AIP Adv. 1, 022103 (2011). 

[25] M. Yamada, T. Hirahara, R. Hobara, S. Hasegawa, 
H. Mizuno, Y. Miyatake, and T. Nagamura, e-J. Surf. 
Sci. Nanotech. 10, 400 (2012). 

[26] X. Wu, X. Li, Z. Song, C. Berger, and W. A. de Heer, 
Phys. Rev. Lett. 98, 136801 (2007). 

[27] F. V. Tikhonenko, D. W. Horsell, R. V. Gorbachev, and 
A. K. Savchenko, Phys. Rev. Lett. 100, 056802 (2008). 

[28] Y.-F. Chen, M.-H. Bae, C. Chialvo, T. Dirks, 
A. Bezryadin, and N. Mason, J. Phys.: Condens. Matter 
22 , 205301 (2010). 

[29] L. Aslamasov and A. Larkin, Phys. Lett. A 26, 238 
(1968). 

[30] M. Tinkham, Introduetion to supereonduetivity (Courier 
Corporation, 2012). 

[31] R. D. Chaudhari and J. B. Brown, Phys. Rev. 139, A1482 
(1965). 

[32] R. Shimizu, K. Sugawara, K. Kanetani, K. Iwaya, 
T. Sato, T. Takahashi, and T. Hitosugi, Phys. Rev. Lett. 
114, 146103 (2015). 



